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Abstract. The frequencies of a cryogenic sapphire oscillator and a hydrogen maser 
are compared to set new constraints on a possible violation of Lorentz invariance. 
We give a detailed description of microwave resonators operating in Whispering 
Gallery modes and then apply it to derive explicit models for Lorentz violating 
effects in our experiment. Models are calculated in the theoretical framework of 
Robertson, Mansouri and Sexl (RMS) and in the standard model extension (SME) 
of Kostelecky and co-workers. We constrain the parameters of the RMS test theory 
to 1/2 - 0ms + S MS = (1.2 ± 2.2) x 10" 9 and [3 MS - a M s - 1 = (1.6 ± 3.0) x 10" 7 
which is of the same order as the best results from other experiments for the former 
and represents a 70 fold improvement for the latter. These results correspond to 
an improvement of our previously published limits [Wolf P. et al., Phys. Rev. Lett. 
90, 6, 060402, (2003)] by about a factor 2. 

1 Introduction 

The Einstein equivalence principle (EEP) is at the heart of special and general 
relativity [1] and a cornerstone of modern physics. One of the constituent 
elements of EEP is Local Lorentz invariance (LLI) which, loosely stated, 
postulates that the outcome of any local test experiment is independent of 
the velocity of the freely falling apparatus (the fundamental hypothesis of 
special relativity) . The central importance of this postulate in modern physics 
has motivated tremendous work to experimentally test LLI [1] . Additionally, 
nearly all unification theories (in particular string theory) violate the EEP 
at some level [2,3] which further motivates experimental searches for such 
violations of the universality of free fall [4] and of Lorentz invariance [5,6]. 

Numerous test theories that allow the modeling and interpretation of ex- 
periments that test LLI have been developed. Kinematical frameworks [7,8] 
postulate a simple parametrisation of the Lorentz transformations with ex- 
periments setting limits on the deviation of those parameters from their spe- 
cial relativistic values. A more fundamental approach is offered by theories 
that parametrise the coupling between gravitational and non-gravitational 
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fields (THe^i [9,1,10] or %g [11] formalisms) which allow the comparison of 
experiments that test different aspects of the EEP. Formalisms based on 
string theory [2-4] have the advantage of being well motivated by theo- 
ries of physics that are at present the only candidates for a unification of 
gravity and the other fundamental forces of nature. Fairly recently a gen- 
eral Lorentz violating extension of the standard model of particle physics 
(Standard Model Extension, SME) has been developed [5] whose Lagrangian 
includes all parametrised Lorentz violating terms that can be formed from 
known fields. Many of the theories mentioned above are included as special 
cases of the SME [12]. 

We report here on experimental tests of LLI using a cryogenic sapphire 
oscillator and a hydrogen maser. After a detailed description of microwave 
resonators operating in Whispering Gallery modes (Sect. 2) we explicitly 
calculate models for Lorentz violating effects in our experiment in the theo- 
retical frameworks of [8] (Sect. 3.1) and the SME (Sect. 3.2). In both cases 
the relative frequency of the two clocks is modulated with, typically, sidereal 
and semi-sidereal periods due to the movement of the lab with the rotation 
of the Earth. The experimental results searching for those variations are pre- 
sented and compared to previously published ones in Sect. 4. We set limits 
on parameters that describe such Lorentz violating effects, improving our 
previous results [13] by a factor 2. Our limits are of the same order as the 
best results from other experiments for the Michclson-Morley type test [15], 
and correspond to a 70-fold improvement over the best other results [14] for 
the Kcnnedy-Thorndike type test. 

2 Whispering Gallery mode resonators 
2.1 Basic principle 

Whispering Gallery (WG) modes in spherical and cylindrical dielectric and 
cavity resonators have been long studied, and have many applications rang- 
ing from novel filter designs to high-Q resonators for oscillator applications 
[19-24]. The boundary conditions ensure that the mode propagates mainly 
around the azimuth in the equatorial plane (in both directions in the case of 
a standing wave). A 'pure' WG mode, as shown in Fig. 1 may be defined with 
only propagation in the azimuthal direction. Thus, for the 'pure' WG mode, 
the Poynting vector has only an component. The mode can be clasified as 
Transverse Electric (TE) or Transverse Magnetic (TM) . TE modes have only 
E r and H z components, while TM modes have only H r and E z components. 

For a dielectric WG mode resonator electromagnetic energy propagates 
just inside the dielectric-free space interface, due to total internal reflections, 
and if we assume the total internal reflection is perfect then we may assume 
100 % of the energy is within the dielectric. Boundary conditions dictate 
that only certain frequencies of resonance occur when an integral number of 
wavelengths fit within the resonator. This happens when 
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Fig. 1. Schematic of the propagation of a 'pure WG' mode in a cylindrical sym- 
metric structure of radius rwG- 



2nr W gy/e t c 

where m is an integer, which represents the number of wavelengths that 
fit along the circumference of the resonator, c is the speed of light in vacuum, 
e the relative permittivity of the dielectric, k m a constant, and t c the travel 
time (in vacuum) of a light signal around the circumference of the resonator. 

This relation is only true for 'pure' WG modes, which only occur as m — > 
oo. At lower values of m the path length is extended by reflections from an 
internal caustic surface (see Fig. 2), and in fact the wave propagates like a 
guided wave around the azimuth. In this case Maxwell's equations must be 
solved, and to calculate the mode frequencies accurately, the roots of the 
corresponding transcendental equation must be found [19]. For this case the 
modes have all electric and magnetic field components, and are classified by 
the dominant field components. Quasi-TE modes have dominant E r and H z 
field components, and are usually denoted as WGE or H z modes. Conversely, 
modes with dominant H r and E z are quasi-TM and are usually denoted as 
WGH or E z modes. Typically, modes are represented with two additional 
numbers n and p, which describe the number of zero crossings (nodes) of 
the dominate field in the radial and axial directions (WG m .„. p ). Thus, the 
fundamental quasi-TE and quasi-TM mode families are written as WGE TO; o,o 
and WGH m; o,o respectively as n = p = 0. 



2.2 Fields in the sapphire resonator 

This subsection describes the techniques to calculate the fields in the uniaxial 
anisotropic sapphire resonator, with the crystal c-axis aligned with the z- 
axis [19]. A schematic of the sapphire resonator is shown in Fig. 3. Two 
techniques were used. The first was a Separation of Variables (SV) technique. 
The advantage was the ability to derive simple analytical expressions for 
all the field components, E z , H z , E^, H$, E r , H r . However, the technique is 
approximate, as it is impossible to calculate consistent analytical expressions 
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Fig. 2. Magnitude of the H z field calculated for the WGEi4 i o,o mode in a sapphire 
disk resonator using the separation of variables technique. The r-cf> and r-z planes 
are represented. The inner caustic is shown, and the mode can be approximated as 
two guided waves propagating in opposite directions around the azimuth. 



for the six field components for an anisotropic right cylinder, with correct 
boundary conditions that satisfy Maxwell's equations. Furthermore, we have 
to make another approximation due to the sloped sides of our resonator (see 
Fig. 3). No doubt the most accurate technique is a numeric method, such 
as Finite Element (FE) analysis. However, for WG modes, field components 
calculated by FE and SV analysis are compared and shown to be close. Also, 
in Sect. 3.2, it is shown that calculations of the sensitivity to Lorentz violation 
in the Standard Model Extension, does not change significantly, whether 
we use the SV or FE results 1%). Even if we use the pure WG mode 
approximation, sensitivity calculations are only over estimated by « 5%. 

Separation of variables The Separation of Variables technique assumes 
the resonator is a perfect cylinder of uniaxial anisotropy with the crystal C- 
axis aligned along the cylinder axis and suspended in free space as shown in 
Fig. 4 [19]. 

The H z magnetic field is symmetric along the z-axis (around z = 0). Given 
these conditions one can show that the solutions to Maxwell's equations yield: 

E z i = iE z ^J m (k E r)cos{m4>)sva.{j3z) 

E z2 = iE Z Q 2 K m (k out r)cos(m4>)sin((3z) 

E z3 = iE z0 3J m (k E r)cos{m4>)e- az (2) 

H z i = H z0 J. m (k H r)sm(m<f>)cos(f3z) 

H z2 = H z02 K m (k out r)sm.{m(j))cos(j3z) 
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Fig. 3. Calculated dimensions of the sapphire cavity after contraction from room 
temperature to 4 K. 




Fig. 4. Solution regions for separation of variables in cylindrical coordinates. 



H z3 = H z03 J m (k H r)sm(m(j))e az 

where the constants and amplitudes are given in tables 1 and 2 below. One 
could also consider a region 4 in the corners where the field has K m (k out r) 
and e~ az dependence. However it is small, and in this work we ignore any 
field there. Also the azimuthal dependence is arbitrary, as long as we keep the 
E and H fields orthogonal either a cos(m^) or sin(m^) dependence could be 
assumed. In actual fact WG modes experimentally exhibit a doublet struc- 
ture due to this degeneracy. In practice the degeneracy is split by a small 
fraction due to imperfections in cylindrical symmetry. Once the z compo- 
nents of the field are calculated, all the other components may be calculated 
from Maxwell's relationships in cylindrical co-ordinates (as shown in refer- 
ence [19]). For example, the fields inside the sapphire (region 1 in Fig. 4) 
which concentrate over 98% of the total energy (c.f. Tab. 4) take the form 



H z = H z0 J m (k H r)sm(m(j))cos(f3z) 
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E z = iE z0 J m (k E r)cos(m(j))sm([3z) 



(3) 



if r = /i rO -^m(fc£;r)sin(m0)sm(/3z) + iJ r0 (J m _i(fc ff r) - J m+ i(fc ff r)) sin(m(/>)sm(/3z) 
r 

_E r = ie r o — J m (fcijr)cos(m0)cos(/3z) + i£yo (■Jm-i(&E7") — Jm+i(^E?")) cos(m0)cos(/3z) 
r 

#0 = h^ -J m {k H r)cos(mcf))su\{(5z) + H^q ( J m -x{k E r) - J m +i(k E r)) cos(m<j))sm(f3z) 
E<p = ie^a-J m {kEr)sva{rn4>)cos{l3z) + iE^o (J m -i(k H r) - J m +i(k H r)) sin(m<£)cos(/32) 

where the amplitudes {E z q, h r0 , H r0 , e r o, E r0 , h^Q, H^q, e^Q, E^q) of the indi- 
vidual terms can be expressed as functions of the constants in Tab. 1 and 
one of the amplitudes (e.g. H z0 ). They are given in Tab. 3 for our resonator. 

An important point must be raised about matching boundary conditions 
to calculate the resonance frequency. The matching at the radial boundary 
generates hybrid solutions where both H z and E z field can co-exist. However, 
it is impossible to simultaneously match all field components on the axial 
boundary, so in some cases approximations must be made. However, this is a 
small effect as only some of the non-dominant components remain unmatched. 

Comparison of separation of variables with finite element analysis 

Due to the sloped sides the SV analysis was checked using FE analysis [25]. 
The calculated frequency was only out by 3 MHz compared to the measured 
frequency of 11.932 GHz. Typically up to 10 MHz discrepancy is normal for 
these types of calculations. The electric field density plot is shown in Fig. 5. 
It is apparent that we should be able to approximate the crystal as a right 
cylinder because the mode is mainly at the perimeter. 




Fig. 5. Electric field density plot inside of the WGEi4 i o,o mode in the r-z plane. 



Given that the permittivity at 4 K is approximately en = 11.349 and 
e_i_ = 9.272, to obtain a frequency close to the measured value and finite 
element calculation, an equivalent height of 21.45 mm is necessary, which 



Whispering Gallery Resonators and Tests of Lorentz Invariance 



7 



gives a frequency of, / = 11.931 GHz close to the measured value of 11.932 
GHz. The solution may be used to calculate the necessary parameters in (2), 
which are summarized in Table 1. 

Table 1. Calculated parameters (in mT 1 ) for the WGEi4 i o i o mode using the SV 
technique. 

ma [3 k H k E k out 

14 707.420 129.571 750.314 830.109 i213.868 



The amplitude of the fields may also be calculated. Because the WG 
mode is dominated by the H z component, we chose to relate the amplitude 
of all other components with respect to H z q. Results are given in Table 2. 
To obtain E z qz uniquely, the component of the tangential E field (E r and 
E$ ) generated by E z , must be assumed to be equal to zero to generate a 
consistent solution (i.e. only assume the H z component exists). 



Table 2. Calculated amplitudes of the field components in (2) with respect to H z o- 





E Z 02 


Ez03 


H Z () 


H Z 02 


Hz03 


19.4538# z0 


6.60905 xlO^Hzo 


428375# z0 


Hzo 


-4.01697xl0~ b ff z0 


355.373# z0 



All other fields arc then obtained from Mawxell relations. For example, 
inside the resonator (region 1 in Fig. 4) they are given by (3) with the am- 
plitudes given in Tab. 3. 



Table 3. Calculated amplitudes of the field components in (3) with respect to H z o- 



Ezo 


h r 


H r o 


e r o 


ErO 


19.4538 H z0 


0.002978 H z0 


-0.08634 Hzo 


-2.344 Hzo 


1.858 Hzo 




H^o 


e<f,o 


E^,o 




-0.003222 H z0 


0.08828 H z0 


-0.06268 Hzo 


62.80 H z0 





Once the amplitude relationships are calculated the percentage of electric 
(Pei) and magnetic field (Prrii) inside and outside the resonator may be 
calculated from: 

Pr = Jv Pm — v - 4) 

' J v e(E*-E)d 3 x 1 J yM (H*-H)d3 x W 
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where the subscript i refers to the component of the field (r, <f> or z). The 
calculated filling factors using SV and FE arc compared in Table 4. 



Table 4. Calculated electric and magnetic filling factors in r, <fi, z directions. 



Method 


Pe z 


Pe r 


Pe 


Pm z 


Prrir 


Pm^ 


FE (in sapphire) 


0.007283 


0.8088 


0.1651 


0.9608 


0.005316 


0.02383 


SV (in sapphire) 


0.006813 


0.8088 


0.1548 


0.9608 


0.005336 


0.02388 


FE (in free space) 


0.0001642 


0.01061 


0.008057 


0.005594 


0.0007796 


0.003710 


SV (in free space) 


0.01139 


0.0106 


0.007682 


0.00555 


0.0007885 


0.003616 



In general the mode is hybrid, and the sloped sides changes slightly the 
amplitude of the non dominant field components (all field components exist, 
although for this mode it is dominantly TE). This was verified by comparing 
the right cylinder to the one with sloped sides using FE analysis. To have 
an analytic solution that is consistent with the FE analysis, the amplitudes 
of some of the fields have been adjusted without changing the form of the 
solutions. This was achieved by leaving the amplitude of the dominant com- 
ponents (H z and E r ) the same, and changing all other components (i?^, H r , 
E z , E^) by a unique amount to make the FE and SV solution of Pe r and 
Pm z the same within the sapphire. After this adjustment, it is clear that the 
solutions are similar. The normalizations factors used are shown in table 5 
and are close to unity. 

Table 5. Amplitude normalization factors of the field components. 



Field comp. 


P[</, 


H r 


E z 




Norm. fact. 


0.82432 


0.73801 


0.569366 


1.00069 



The biggest discrepancy is with the Pe z component in free space. This is 
not surprising as it was the E z q% amplitude in (2), which was compromised 
with the SV technique. However, the technique allows us to obtain simple 
analytic solutions of the form (3) to the field for a WG mode resonator that 
can be useful when calculating the sensitivity to Lorcntz violations within 
the framework of the Standard Model Extension [12] (see Sect. 3.2). 

3 Theoretical analysis 

Owing to their simplicity the kinematical frameworks of [7,8] have been 
widely used to model and interpret many previous experiments testing LLI 
[13,14,16,18,26,27] and we will follow that route in Sect. 3.1. An analysis 
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based on the more fundamental "Standard Model Extension" (SME) [5,12] 
is then presented in Sect. 3.2. 

By construction, kinematical frameworks do not allow for any dynamical 
effects on the measurement apparatus. This implies that in all inertial frames 
two clocks of different nature (e.g. based on different atomic species) run at 
the same relative rate, and two length standards made of different materials 
keep their relative lengths. Coordinates are defined by the clocks and length 
standards, and only the transformations between those coordinate systems 
are modified. In general this leads to observable effects on light propagation 
in moving frames but, by definition, to no observable effects on clocks and 
length standards. In particular, no attempt is made at explaining the un- 
derlying physics (e.g. modified Maxwell and/or Dirac equations) that could 
lead to Lorentz violating light propagation but leave e.g. atomic energy levels 
unchanged. On the other hand dynamical frameworks (e.g. the THe/x formal- 
ism or the SME) in general use a modified general Lagrangian that leads to 
modified Maxwell and Dirac equations and hence to Lorentz violating light 
propagation and atomic properties, which is why they are considered more 
fundamental and more complete than the kinematical frameworks. Further- 
more, as shown in [12], the SME is kept sufficiently general to, in fact, en- 
compass the kinematical frameworks and some other dynamical frameworks 
(in particular the THe/i formalism) as special cases, although there are no 
simple and direct relationships between the respective parameters. 

3.1 The Robertson, Mansouri &: Sexl framework 

Kinematical frameworks for the description of Lorentz violation have been 
pioneered by Robertson [7] and further refined by Mansouri and Sexl [8] 
and others. Fundamentally the different versions of these frameworks are 
equivalent, and relations between their parameters are readily obtained. As 
mentioned above these frameworks postulate generalized transformations be- 
tween a preferred frame candidate E(T, X) and a moving frame S(t, x) where 
it is assumed that in both frames coordinates are realized by identical stan- 
dards. We start from the transformations of [8] (in differential form) for the 
case where the velocity of 5* as measured in £ is along the positive X-axis, 
and assuming Einstein synchronization in S (we will be concerned with signal 
travel times around closed loops so the choice of synchronization convention 
can play no role): 

, m 1 / , vdx\ dx v ( , vdx\ „„ dy dz 

dT=-(dt + —y,dx = - + -(dt+—y,dY=^dz=- (5 ) 

with c the velocity of light in vacuum in S. Using the usual expansion of the 
three parameters (a w 1 + a M s« 2 /c 2 + 0(4); b w 1 + /3 M s« 2 /c 2 + 0(4); d w 
l+4isv 2 /c 2 +0(4)), setting c 2 dT 2 = dX 2 +dY 2 +dZ 2 in E, and transforming 
according to (5) we find the coordinate travel time of a light signal in S: 
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dt=^(l- (/3 MS - a M s ~ 1) ^ - (\ ~ ^ MS + Sms ) sin2<? ^2 ) + °^ 

(6) 

where dl = \J dx 2 + dy 2 + dz 2 and 9 is the angle between the direction of 
light propagation and the velocity v of S in S. In special relativity qms = 
— 1/2; /3ms = 1/2; 5ms = and (5) reduces to the usual Lorcntz transforma- 
tions. Generally, the best candidate for S is taken to be the frame of the 
cosmic microwave background (CMB) [28,29] with the velocity of the solar 
system in that frame taken as Vq w 377 km/s, decl. w —6.4°, RA w 11. 2h. 

Michelson-Morley type experiments [30,16] determine the coefficient Pmm — 
(1/2 — /3ms + <$ms) of the direction dependent term. For many years the most 
stringent limit on that parameter was |Pmm| < 5 x 10~ 9 determined over 
23 years ago in an outstanding experiment [16]. Our experiment confirms 
that result with roughly equivalent uncertainty (2.2 x 10~ 9 ). Recently an im- 
provement to | Pmm | < l-5x 10~ 9 has been reported [15]. Kennedy-Thorndikc 
experiments [17,18,14] measure the coefficient Pkt — (/3ms ~ «ms ~ 1) of the 
velocity dependent term. The most stringent limit [14] on \Pkt\ has been 
recently improved from [18] by a factor 3 to \Pkt\ < 2.1 x 1CT 5 . We improve 
this result by a factor of 70 to |Pr-t| < 3.0 x 10~ 7 . Finally clock comparison 
and Doppler experiments [26,27,31] measure «ms- The most stringent result 
comes from the recent experiment of [32], limiting |o:ms + 1/2| to < 2.2x 10 -7 . 
The three types of experiments taken together then completely characterize 
any deviation from Lorentz invariance in this particular test theory, with 
present limits summarized in Table 6. 



Table 6. Present limits on Lorentz violating parameters in the framework of [8] 



Reference a M s + 1/2 1/2 - /? M s + 5ms (3ms - «ms - l 

[26,31,27] < 8 x 1(T 7 

[32] < 2.2 x 1(T 7 

[16] - < 5 x 1(T 9 

[15] - (2.2 ± 1.5) x 1(T 9 

[14] - - (1.9 ±2.1) x 1(T 5 

our previous results [13] - (1.5 ± 4.2) x 1(T 9 (-3.1 ± 6.9) x 1(T 7 

this work - (1.2 ± 2.2) x 1(T 9 (1.6 ± 3.0) x 10~ 7 



As described in Sect. 2 our cryogenic oscillator consists of a sapphire 
crystal of cylindrical shape operating in a whispering gallery mode, and its 
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coordinate frequency can be expressed by equation (1) where t c is the coor- 
dinate travel time of a light signal around the circumference of the cylinder. 
Calculating t c from (6) the relative frequency difference between the sapphire 
oscillator and the hydrogen maser (which, by definition, realizes coordinate 
time in S [33]) is 

+ PMMfi r S in^,^ + C>(3) (7) 
i/q c c In J 

where v$ is the unperturbed frequency, v(t) is the (time dependent) speed of 
the lab in S, and ip is the azimuthal angle of the light signal in the plane of 
the cylinder. The periodic time dependence of v and 9 due to the rotation 
and orbital motion of the Earth with respect to the CMB frame allow us to 
set limits on the two parameters in (7) by adjusting the periodic terms of 
appropriate frequency and phase (see [35] for calculations of similar effects 
for several types of oscillator modes) . Given the limited durations of our data 
sets (< 15 days) the dominant periodic terms arise from the Earth's rotation, 
so retaining only those we have v(t) = u + u> x R with u the velocity of the 
solar system with respect to the CMB, ui the angular velocity of the Earth, 
and R the geocentric position of the lab. We then find after some calculation. 

Av/vq = Pkt{Hsvo.X) 

+Pmm{Acos\ + i?cos(2A) + CsinA + DsinAcosA + £sinAcos(2A)) 

(8) 

where A = u>t + <fi, and A-E and <f> are constants depending on the latitude 
and longitude of the lab (w 48.7°N and 2.33°E for Paris). Numerically H w 
-2.6 x 1CT 9 , A w -8.8 x 1(T 8 , B w 1.8 x 1CT 7 , C-E of order 1(T 9 . We 
note that in (8) the dominant time variations of the two combinations of 
parameters are in quadrature and at twice the frequency which indicates 
that they should decorelate well in the data analysis allowing a simultaneous 
determination of the two (as confirmed by the correlation coefficients given 
in Sect. 4). Adjusting this simplified model to our data we obtain results 
that differ by less than 10% from the results presented in Sect. 4 that were 
obtained using the complete model ((7) including the orbital motion of the 
Earth). 

3.2 The Standard Model Extension 

The fundamental theory of the Standard Model Extension (SME) as applied 
to electrodynamics is laid out in [12]. Here we summarise the main points 
relating to e-m cavity tests, and apply them to model our experiment. 

The photon sector of the SME is described by a Lagrangian that takes 
the form 
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C = --F^ v + \{k AF ) K e K ^ v A x F^ - \{k F ) K ^ v F KX F^ (9) 

where F M „ = d^A^ — d v A^. The first term is the usual Maxwell part while the 
second and third represent Lorentz violating contributions that depend on the 
parameters Uaf and kp. For most analysis the Uaf parameter is set to for 
theoretical reasons (c.f. [12]) but which is also well supported experimentally. 
The remaining dimensionless tensor (kp)^^ has a total of 19 independent 
components that need to be determined by experiment. Retaining only this 
term leads to Maxwell equations that take the familiar form but with D and 
H fields in vacuum defined by a six dimensional matrix equation 



D\ / £ o(l + kde) 



kdb 



(1 + «j?b) 




(10) 



where the k are 3x3 matrices whose components are particular combinations 
of the Uf tensor (c.f. equation (5) of [12]). Equation (10) indicates a useful 
analogy between the SME in vacuum and standard Maxwell equations in 
homogeneous anisotropic media. For the analysis of different experiments it 
turns out to be useful to introduce further combinations of the k matrices 
defined by: 







^{k de + K HB y k 


(«e- 


_ 


t^{kde - K HB y k 


(a 0+ . 


= 


-{k db + n HE ) jk 


(«o- 


pfc = 


^{kdb - n HE y k 


{Rtr. 


)jk _ 





(11) 

The first four of these equations define traceless 3x3 matrices, while 
the last defines a single coefficient. All k matrices are symmetric except k Q+ 
which is antisymmetric. These characteristics leave a total of 19 independent 
coefficients of the k. In general experimental results are quoted and compared 
using the k parameters rather than the original ftp tensor components and 
this is the route we will follow in the present analysis. 

The kp tensor in (9), and consequently the k tensors in (10) and (11), are 
frame dependent and consequently vary as a function of the coordinate system 
chosen to analyse a given experiment. In principle they may be constant and 
non-zero in any frame (e.g. the cavity frame or the lab frame). However, 
any non-zero values are expected to arise from Planck-scale effects in the 
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early Universe. Therefore the components of fc^ should be constant in a 
cosmological frame (e.g. the one defined by the CMB radiation) or any frame 
that moves with a constant velocity and shows no rotation with respect to the 
cosmological one. Consequently the conventionally chosen frame to analyse 
and compare experiments in the SME is a sun-centred, non-rotating frame 
as defined in [12]. The general proceedure is to calculate the perturbation of 
the resonator frequency as a function of the unperturbed E n and B fields 
and k tensors in the lab frame and then to transform the n tensors to the 
conventional sun-centred frame. This transformation will introduce a time 
variation of the frequency related to the movement of the lab with respect to 
the sun-centred frame (typically introducing time variations of sidereal and 
semi-sidereal periods for an Earth fixed experiment). 

In [12] the authors derive an expression for the perturbed frequency of a 
resonator of the form 



— = - J- / d 3 x (e E * • k de • E - Mo 1b o* ■ «hb ■ B (12) 
v o \U) Jv 

+ 2Rc( A /^E * -KDB -B )) 
V Mo / 

where Bo,H ,Eo,D are the unperturbed (standard Maxwell) fields and 
(U) = J v d 3 x(Eo ■ Do* + Bo • Ho*). This expression can be applied directly 
to our resonator using the fields calculated in Sect. 2.2. 

The resonator is placed in the lab with its symmetry axis along the verti- 
cal. So applying (13) in the lab frame (z-axis vertical upwards, x-axis pointing 
south), with the fields calculated as described in Sect. 2.2, we obtain an ex- 
pression for the frequency variation of the resonator 



+ (M H B)t: b ((KH B )t: b + {k hb )%) + (M HB )i: b (K HB )i: b as) 

with the M.i ab components given in Tab. 7. To obtain the values in Tab. 7 
we take into account the fields inside the resonator (e.g. (3)) and outside 
(< 2% of the energy). We note that when using coordinates in which the 
K \a b components are constants over the volume of the cylinder (e.g. cartesian 
coordinates), the integrals in (13) are equivalent to filling factor integrals 
(c.f. (4)) for the diagonal terms. This simplifies the calculation of the Adf ab 
components when using numerical techniques, as simple algorithims exist to 
calculate filling factors without explicitly calculating the field [36]. 

The last step is to transform the k tensors in (13) to the conventional 
sun-centred frame using the explicit transformations provided in [12], and to 
express the result in terms of the k tensors of (11). We obtain 
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Table 7. Mi a b components calculated using the different techniques for the deter- 
mination of the fields inside the resonator described in Sect. 2.2. 

M]j b FE SV PurcWG 

(MDE)fab -0.03093 -0.0355 -0.02696 

(M D E) z iab -0.0004030 -0.005996 0.00 

(MHB)fab 0.008408 0.008405 0.00 
{M HB )T a b 0-4832 0.4832 0.50 



— = K + C A cos(QT) + S A sin(i?T) 
+ C w cos(wT e ) + S w sin(o;T e ) 

+ C 2a) cos(2^T ffi ) + 5 2w sin(2wT e ) (14) 
+ C w+ r2Cos((w + J?)T© + tp) + S^+nsin^uj + J?)T© + tp) 
+ C u -ncos((w - Q)T & — <p) + S u ,_^sin(((o; - J?)T e - <p) 
+ C 2w+f2 cos((2w + /2)T e + tp) + S 2u +nsm((2u + tt)T ffi + <p) 
+ C 2w _ r2 cos((2w - (2)T e -<p) + 5 , 2w _r 2 sin((2w - (2)T e - tp) 

where u and Q are the sidereal and anual angular frequencies (w = 7.292115 x 
10" 5 rad/s, = 1.991 x 10" 7 rad/s), T e is the time since a coincidence of 
the lab y axis with the sun-centred Y axis (as defined in [12]), T is the time 
since a spring equinox and ip is a constant given by tp = S7(T — T©). 

The K e + and k - tensors have been determined to < 2 x 1CP 32 by astro- 
physical tests [12], orders of magnitude below what we can hope to achieve 
in resonator experiments. Consequently we set those parameters to zero, ob- 
taining the C and S coefficients of (14) given in Tab. 8. They involve 4 com- 
binations of the 5 independent components of K e _ and all three independent 
components of k Q+ , but do not involve k tr . So our experiment is sensitive to 
7 of the 9 parameters not determined by astrophysical tests. The numerical 
values in Tab. 8 were obtained using the FE results in Tab. 7. Similar values 
are obtained when using the SV or 'pure WG' results, with the differences 
not exceeding 1% and 5% respectively. This small discrepancy is due to the 
fact that the dominant term ((.Mhb)^) m Tab. 7 is almost unaffected by 
the technique used. However, after transformation to the sun centred frame 
and to the k tensors it is that term which dominates the sensitivity to the k 
components. 

Present limits on the 7 parameters are summarised in Tab. 9. At least 
one year of regular data is required to be able to decorrelate the ui and 2lo 
frequencies from the to ± fi and 2oj ± Q ones. The data of our experiment 
that we present here (see next section) are not yet sufficient for that purpose 
(correlation coefficients are still w 0.4), so the next section only presents the 
results in the kinematical framework of Sect. 3.1. However, the experiment is 
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still running, and we expect to publish our results on the SME test by early 
2004. 

Table 8. Coefficients in (14). The constant K and the annual coefficients (Ca and 
Sa) are not given as our experiment is insensitive to those terms. 



frequency C S 

uj -0.4402 kf_ + (1.066 x 10~ b ) k^ + -0.4402 k'f_ ± (1.066 x 10~ b ) kf + 

2uj -0.09625 (ft" - Hf-) -0.1925 kf_ 

cu + n -(8.6405 x 10" 6 ) kf + (8.6405 x 10~ 6 ) kf + + (1.779 x 10~ 6 ) kJ ; 

w - O -(8.6405 x 10~ 6 ) kf + (8.6405 x 10~ 6 ) kf + - (4.196 x 10 -5 ) k^ 

2uj + n (7.780 x 10~ 7 ) kf + (7.780 x 10~ 7 ) kf+ 

2lu - Q -(1.834 x IP" 5 ) kjl -(1.834 x 10~ 5 ) kf + 



Table 9. Present limits (la uncertainties) on the 7 SME parameters our experiment 
is sensitive to, as determined in [15]. Values are given in 10~ 15 for k e - and 10~ n 
for k a+ . 



K e '_ 




~±i -li ~-±-± 








1.7 ±2.6 


-6.3 ± 12.4 


3.6 ±9.0 8.9 ±4.9 


-14 ± 14 


1.2 ±2.6 


-0.1 ±2.7 



We conclude this section by a short discussion of the effects of non-zero 
values of the SME parameters on our frequency reference (the hydrogen 
mascr) and the crystal structure of the sapphire (and hence the resonator 
shape). Such effects would add to the direct effect on the electromagnetic 
fields and may therefore lead to overall cancellation or enhancement of the 
sensitivities calculated above. However, it was shown in [37] that the effect 
on the sapphire crystal amounts to only a few percent of the direct effect 
on the fields, and [38] show that the hydrogen mp = — > m' F = clock 
transition is not affected to first order. Hence the total effect is dominated 
by the Lorentz violating properties of the electromagnetic fields inside the 
resonator and well described (to a few percent) by the model derived above. 



4 Experimental results 

As mentioned in the previous section, the experimental data presented here 
does not yet allow a complete decorrelation of the different parameters in the 
SME model, so we concentrate in this section on the experimental results in 
the kinematic (RMS) framework presented in section 3.1. 
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The cryogenic Sapphire Oscillator (CSO) is an active system oscillating at 
the resonant frequency (i.e. a classical loop oscillator which amplifies and re- 
injects the "natural" resonator signal). Additionally the signal is locked to the 
resonance using the Pound-Drever technique (modulation at w 80 kHz). The 
incident power is stabilised in the cryogenic environment and the spurious 
AM modulation is minimised using a servo loop. To minimise temperature 
sensitivity the resonator is heated (inside the 4 K environment) and stabilised 
to the temperature turning point (« 6 K) of the resonator frequency which 
arises due to paramagnetic impurities in the sapphire. Under these conditions 
the loaded quality factor of the resonator is slightly below 10 9 . The resonator 
is kept permanently at cryogenic temperatures, with helium refills taking 
place about every 20 - 25 days. 

The CSO is compared to a commercial (Datum Inc.) active hydrogen 
maser whose frequency is also regularly compared to caesium and rubidium 
atomic fountain clocks in the laboratory [39]. The CSO resonant frequency at 
11.932 GHz is compared to the 100 MHz output of the hydrogen maser. The 
maser signal is multiplied up to 12 GHz of which the CSO signal is subtracted. 
The remaining w 67 MHz signal is mixed to a synthesizer signal at the same 
frequency and the low frequency beat at w 64 Hz is counted, giving access to 
the frequency difference between the maser and the CSO. The instability of 
the comparison chain has been measured at < 2 x 10 _14 r _1 , with long term 
instabilities dominated by temperature variations, but not exceeding 10~ 16 . 

Since September 2002 we are taking continuous temperature measure- 
ments on top of the CSO dewar and behind the electronics rack. Starting 
January 2003 we have implemented an active temperature control of the 
CSO room and changed some of the electronics. As a result the diurnal and 
semi-diurnal temperature variations during measurement runs (« 2 weeks) 
were greatly reduced to less than 0.025° C in amplitude (best case), and 
longer and more reliable data sets became available. 

Our previously published results [13] are based on data sets taken be- 
tween Nov. 2001 and Sep. 2002. Most of that data (except the last data set) 
was taken before implementation of permanent temperature control. As a 
result the uncertainties in [13] were dominated by the systematic effects from 
temperature variations. Here we use only data that was permanently tem- 
perature controlled, 13 data sets in total spanning Sept. 2002 to Aug. 2003, 
of differing lengths (5 to 16 days, 140 days in total), thereby significantly 
reducing the uncertainties from systematic effects. The sampling time for all 
data sets was 100 s except two data sets with To = 12 s. To make the data 
more manageable we first average all points to r = 2500 s. For the data 
analysis we simultaneously adjust an offset and a rate (natural frequency 
drift, typically w 1.7 x 10~ 18 s -1 ) per data set and the two parameters of the 
model (7). In the model (7) we take into account the rotation of the Earth 
and the Earth's orbital motion, the latter contributing little as any constant 
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or linear terms over the durations of the individual data sets are absorbed 
by the adjusted offsets and rates. 

When carrying out an ordinary least squares (OLS) adjustment we note 
that the residuals have a significantly non-white behavior. The power spectral 
density (PSD) of the residuals when fitted with a power law model of the form 
Sy(f) — kf^ yields typically /z w —1.5 to —2. In the presence of non-white 
noise OLS is not the optimal regression method [40,41] as it can lead to 
significant underestimation of the parameter uncertainties [40] . 

An alternative method is weighted least squares (WLS) [41] which al- 
lows one to account for non-random noise processes in the original data by 
pre-multiplying both sides of the design equation (our equation (7) plus the 
offsets and rates) by a weighting matrix containing off diagonal elements. To 
determine these off diagonal terms we first carry out OLS and adjust the 
S y (f) — kf 1 * power law model to the PSD of the post-fit residuals determin- 
ing a value of /x for each data set. We then use these xx values to construct 
a weighting matrix following the method of fractional differencing described, 
for example, in [40]. Figure 6 shows the resulting values of the two parameters 
(Pkt and Pmm) for each individual data set. A global WLS fit of the two 
parameters and the 13 offsets and drifts yields Pmm = (1.2 ± 1.9) x 10^ 9 and 
Pkt — (1-6 ± 2.3) x 10~ 7 (la uncertainties), with the correlation coefficient 
between the two parameters less than 0.01 and all other correlation coeffi- 
cients < 0.06. The distribution of the 13 individual values around the ones 
obtained from the global fit is well compatible with a normal distribution (x 2 
= 10.7 and x 2 = 14-6 for Pmm and Pkt respectively). 

Systematic effects at diurnal or semi-diurnal frequencies with the appro- 
priate phase could mask a putative sidereal signal. The statistical uncertain- 
ties of Pm m and Pkt obtained from the WLS fit above correspond to sidereal 
and semi-sidereal terms (from (8)) of w 7 x 10~ 16 and «4x 10~ 16 respec- 
tively so any systematic effects exceeding these limits need to be taken into 
account in the final uncertainty. We expect the main contributions to such 
effects to arise from temperature, pressure and magnetic field variations that 
would affect the hydrogen maser, the CSO and the associated electronics, 
and from tilt variations of the CSO which are known to affect its frequency. 

To estimate the tilt sensitivity we have intentionally tilted the oscillator 
by w 5 mrad off its average position which led to relative frequency variations 
of s=a 3 x 10 -13 from which we deduce a tilt sensitivity of w 6 x 10 _17 /xrad _1 . 
Measured tilt variations in the lab at diurnal and semi-diurnal periods were 
measured at 4.6 /xrad and 1.6 iirad respectively which leads to frequency 
variations that do not exceed 3 x 10~ 16 and 1 x 10~ 16 respectively and are 
therefore negligible with respect to the statistical uncertainties. 

In December 2002 we implemented an active temperature stabilization 
inside an isolated volume (m 15m 3 ) that included the CSO and all the asso- 
ciated electronics. The temperature was measured continously in two fixed 
locations (behind the electronics rack and on top of the dewar). For the best 
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data sets the measured temperature variations did not exceed 0.02/0.01 °C 
in amplitude for the diurnal and semi-diurnal components. In the worst cases 
(the two 2002 data sets and some data sets taken during a partial air con- 
ditioning failure) the measured temperature variations could reach 0.26/0.08 
°C. When intentionally heating and cooling the CSO lab by sw 3°C we see 
frequency variations of « 5 x 10~ 15 per °C. This is also confirmed when we 
induce a large sinusoidal temperature variation (« 1.5°C amplitude). Using 
this we can calculate a value for temperature induced frequency variations at 
diurnal and semi-diurnal frequencies for each data set, obtaining values that 
range from w 5 x 10~ 17 to w 1.3 x 10~ 15 . 

The hydrogen mascr is kept in a dedicated, environmentally controlled 
clock room. Measurements of magnetic field, temperature and atmospheric 
pressure in that room and the maser sensitivities as specified by the manu- 
facturer allow us to exclude any systematic effects on the maser frequency 
that would exceed the statistical uncertainties above and the systematic un- 
certainties from temperature variations in the CSO lab. 

Our final uncertainties (the error bars in Fig. 6) are the quadratic sums of 
the statistical uncertainties from the WLS adjustment for each data set and 
the systematic uncertainties calculated for each data set from (8) and the 
measured temperature variations. For the global adjustment we average the 
systematic uncertainties from the individual data sets obtaining ±1.2 x 10~ 9 
on Pmm and ±1.9 x 10~ 7 on Pkt- Adding these quadratically to the WLS 
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Fig. 6. Values of the two parameters (Pkt and Pmm) from a fit to each individual 
data set (blue diamonds) and a global fit to all the data (red squares). For compari- 
son our previously published results [13] are also shown (green triangles). The error 
bars indicate the combined uncertainties from statistics and systematic effects. 
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statistical uncertainties of the global adjustment we obtain as our final result 
Pmm = (1-2 ± 2.2) x f(T 9 and P KT = (1.6 ± 3.0) x 10~ 7 (la uncertainties). 

5 Conclusion 

We have presented a detailed description of whispering gallery modes in 
cylindrical resonators. This includes analytical and numerical calculations 
of electro-magnetic fields and energy filling factors for the sapphire resonator 
used in our experiment. It was shown that the different calculations gave 
similar results (differences at the % level). We then applied those results 
to model Lorentz violating frequency shifts of the resonator in two different 
theoretical frameworks: the kinematical framework of Robertson, Mansouri 

6 Sexl (RMS), and the standard model extension (SME) developed by Kost- 
clccky and co-workers. In both cases we obtain explicit expressions ((7) and 
(14)) for the fractional frequency variation of our resonator as a function of 
the Earth's rotational and orbital angular velocities and of the parameters 
describing Lorentz violation in RMS or the SME. We show that the exper- 
imental sensitivity to the SME parameters is not affected by more than 1 
% when using the analytical or numerical method to calculate the resonator 
fields. 

Experimental results were given for the RMS test, showing that our ex- 
periment simultaneously constrains two combinations of the three param- 
eters of the Mansouri and Sexl test theory (previously measured individu- 
ally by Michelson-Morley and Kennedy-Thorndikc experiments). We obtain 
8ms — Pms + 1/2 = 1.2(1.9)(1.2) x 10~ 9 which is of the same order as the 
best previous results [15,16], and /3ms — oms — 1 = 1-6(2. 3)(1. 9) x 1CP 7 which 
improves the best previous limit [14] by a factor of 70 (the first bracket indi- 
cates the la uncertainty from statistics the second from systematic effects). 
We improve our own previous results [13] by about a factor 2 due to more 
and longer data sets and to improved temperature control of the experiment 
(see Tab. 6 for a summary of present limits). We note that our value on 
8ms — /?ms + 1/2 is compatible with the slightly significant recent result of 
[15] who obtained 5 MS - f3 M s + 1/2 = (2.2 ± 1.5) x 10~ 9 . 

As a result of our experiment the Lorentz transformations are confirmed 
in the RMS framework with an overall uncertainty of < 3 x 10~ 7 limited by 
our determination of /3ms — a MS — 1 and the recent limit [32] of 2.2 x 10~ 7 
on the determination of qms- The latter is likely to improve in the coming 
years by experiments such as ACES (Atomic Clock Ensemble in Space [42]) 
that will compare ground clocks to clocks on the international space station 
aiming at a 10 -8 measurement of qms- 

Concerning the SME we do not present experimental results as our data 
is not yet sufficient to decorrelate all SME parameters. However, the exper- 
iment is still running and we expect to have enough data by early 2004 to 
obtain unambiguous valus of all SME parameters. Nonetheless, we can al- 
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ready note that our experiment is complementary with respect to the best 
previous determination of those parameters [15] by looking at tables 8 and 
9. From Tab. 9 we see that the most accurate values are those for h^L and 
— Kg'l , the two parameters our experiment is less sensitive to (c.f. second 
line of Tab. 8) whereas our best sensitivity is to k],_ and (first line of Tab. 
8) which are less well determined in [15]. Given our uncertainties at sidereal 
and semi-sidereal periods of w 1 x 10~ 15 and w 5 x 1CT 16 respectively (c.f. 
Sect. 4) we expect to obtain uncertainties of about 2 x 10~ 15 on k^_, kt^L 
and about 3 — 5 x 10~ 15 on k^_, n lx _ — k^_. Those results, if achieved, will 
compare favourably to the present limits of Tab. 9. 

For the future we do not expect significant improvements using our present 
experimental setup, due to the already relatively long total data span and 
to expected systematic limits from both, the hydrogen maser and the sap- 
phire resonator in the low 1CP 16 region. Significant improvements in the near 
future are more likely to come from new proposals, for example, using two or- 
thogonal resonators or two orthogonal modes in the same sapphire resonator 
placed on a rotating platform [35]. Such a set-up is likely to improve the 
tests of LLI by several orders of magnitude as the relevant time variations 
will now be at the rotation frequency (« 0.01 — 0.1 Hz) which is the range 
in which such resonators are the most stable (~ 100 fold better stability). 
Additionally many systematic effects are likely to cancel between the two 
orthogonal oscillators or modes and the remaining ones are likely to be less 
coupled to the rotation frequency than to the sidereal frequencies used in 
our experiment. Ultimately, it has been proposed [43] to conduct these tests 
on board an Earth orbiting satellite, again with a potential gain of several 
orders of magnitudes over current limits. 
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